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ABSTRACT 

The distortion of flames in flows with vortical motion is examined via asymptotic analysis and 
numerical simulation. The model consists of a constant-density, one-step, irreversible Arrhenius 
reaction between initially unmixed species occupying adjacent half-planes which are then allowed 
to mix and react in the presence of a vortex. The evolution in time of the temperature and mass- 
fraction fields is followed. Emphasis is placed on the ignition time and location as a function of 
vortex Reynolds number and initial temperature differences of the reacting species. The study 
brings out the influence of the vortex on the chemical reaction. In all phases, good agreement is 
observed between asymptotic analysis and the full numerical solution of the model equations. 


The second and third authors were supported by the National Aeronautics and Space Administration under NASA Contract No. NAS1- 
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1. INTRODUCTION 


In a seminal paper, Marble (1985) developed a simple model problem which describes the distortion 
of a diffusion flame by the presence of a viscous vortex flow, as might be found in the large scale structures 
of turbulent diffusion flames. Assuming fast chemical kinetics and the constant density approximation, a 
flame sheet exists at time zero separating a region of fuel in the upper half plane and a region of oxidizer in 
the lower half plane. As time increases, the flame sheet is allowed to evolve by diffusion and convection 
and is wound up by the vortex flow field. The analysis was based on the assumption that the flame elements 
are locally strained and sheared, but that the shearing motion was unimportant and hence could be 
neglected. Theoretical analysis shows that (i) the vortex field increases the reactant consumption, (ii) the 
reaction rate is independent of time, (iii) a similarity rule for the core radius growth exists, and (iv) a simi- 
larity rule exists for the reactant consumption rate. 

The analytical problem posed by Marble (1985) (henceforth called the Marble problem) has been 
solved numerically by Laverdant and Candel (1988) and Rehm, Baum, Lozier and Aronson (1989), but 
without the strain-shear approximation. Thus the complete flow field was considered. Laverdant and Candel 
(1988) solved the Marble problem by an implicit finite difference scheme. They verified numerically 
Marble’s similarity rule for the core radius growth and the result that the reaction rate is independent of 
time, and also verified the correctness of the shear approximation. In addition, the effect of varying the 
stoichiometric (or equivalence) parameter was also studied. They showed that the core is mostly surrounded 
by the oxidizer when the stoichiometric parameter is less than unity (fuel lean), and is mosdy surrounded by 
the fuel when it is greater than unity (fuel rich). Rehm, et al. (1989) solved the Marble problem by first 
transforming the system into a Lagrangian coordinate system. Then assuming a similarity solution, the 
resulting problem was solved analytically using Fourier series. The Fourier amplitudes were determined by 
either a two-point boundary-value solver, valid for all values of the independent parameters, or in the 
asymptotic limit of large Schmidt number (as might be more appropriate for liquids than for gases). The 
results of this study confirmed the behaviour of the numerical solution of Laverdant and Candel (1988), as 
well as the dependency of the reactant consumption rate found by Marble (1986). Norton (1983) also 
solved the Marble problem with the strain-shear approximation, but using finite rate chemistry instead of the 
flame sheet approximation, confirming the results of Marble. 

Many extensions of the Marble problem exist, each utilizing the same analytical approach as that of 
Marble (1985) and assuming the fast chemistry limit, some of which are briefly reviewed here. Karagozian 
and Marble (1986) considered the time dependent interaction of a diffusion flame with an axially strained 
vortex in the third dimension. It was shown that the main conclusions for the three dimensional vortex field 
remains formally identical to the two dimensional vortex field (the Marble problem). In addition, density 
effects were shown to be small. Karagozian and Manda (1986) considered the distortion of a two dimen- 
sional fuel strip in the presence of a pair of counter-rotating vorticies. The fuel strip was assumed either 
infinite or semi-infinite in extent, separated from the oxidizer by two flame sheets. On each flame sheet 
resides a vortex, which are allowed to rotate in opposite directions. The semi-infinite fuel strip is particu- 
larly relevant to the vortical flame structure formed at a circular orifice or nozzle. Cetegen and Sirignano 
considered the interaction of a diffusion flame with a single vortex (1987) and the interaction with an 
infinite row of two dimensional vortices, which is representative of the temporal growth of a mixing layer 
before pairing has occurred (1988). In both cases, the emphasis was on the construction of the concentration 
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profilcs in the vortical structure, as well as, the construction of the probability density functions (pdfs). A 
review of these two cases can be found in Sirignano (1989). Peters and Williams (1988) developed an 
analogous problem for a premixed flame, relevant for describing premixed turbulent combustion in large 
scale structures. The emphasis of this study was in describing the growth of the reacted core, flame extinc- 
tion by stretch for Lewis numbers greater than unity, and the effect of heat release. 

In another seminal paper, Linan and Crespo (1976) examined the structure of a diffusion flame in the 
unsteady mixing of two half spaces of fuel and oxidizer. By using a combination of large activation energy 
asymptotics and numerics, they analyzed the continuous evolution from nearly frozen flow to near equili- 
brium flow and showed that three laminar regimes exist within the flow; the ignition, deflagration, and 
diffusion flame regimes. The ignition regime is a region where the combustible gases mix until, at some 
finite time, a thermal explosion occurs at a well defined location and the gas is ignited. The second regime 
is the deflagration region. After ignition, a pair of deflagration waves (or premixed flamelets ) emerge 
according to classical thermal explosion theory. These waves arise because the mixture is pot stoichiometric 
in the premixed region. One of the reactants is consumed locally, leaving behind an excess of the other 
reactant. Thus, one of the flamelets is fuel-rich and the other is fuel-lean. There is excess fuel behind the 
fuel-rich flamelet and excess oxidizer behind the fuel-lean flamelet. Concentration gradients behind the 
flamelets drive the unbumt fuel and oxidizer towards the diffusion flame where they are consumed locally. 
These premixed flamelets are quite weak in that the temperature rise associated with them is small, and they 
exist only until all of the deficient reactant is consumed. Just beyond the deflagration waves, a diffusion 
flame regime exists where the mixing process is governed by diffusion in the direction normal to the flame. 
As time increases, the diffusion flame approaches a flame sheet. We note that the existence of a well 
defined ignition point and the premixed flamelets depends critically on the relative magnitudes of the two 
initial temperatures to that of the adiabatic diffusion flame temperature. That is, if the adiabatic diffusion 
flame temperature is greater than either initial temperature, a well defined ignition point always occurs, fol- 
lowed by the premixed flamelets. On the other hand, if the adiabatic diffusion flame temperature is between 
the two initial temperatures, there is no well defined ignition point, and a single premixed flame merges 
smoothly into the diffusion flame. The analytical results of Linan and Crespo were extended by Jackson and 
Hussaini (1988) and Grosch and Jackson (1991) in their studies of steady supersonic combustion in 
compressible laminar mixing layers. In particular, the three flame regimes also exist at high speeds provided 
the adiabatic diffusion flame temperature is greater than either freestream temperature. Finally, we note that 
this configuration of two flamelets and a diffusion flame (collectively called a tribrachial flame) has also 
been observed in the propagation of a premixed flame into a nonuniform mixture (e.g., Buckmaster and 
Matalon, 1988). 

Our goal is to combine the models of Marble (1985) and Linan and Crespo (1976). That is, we pro- 
pose to investigate the continuous evolution of initially unmixed species occupying adjacent half planes to 
near equilibrium flow in the presence of a viscous vortex without invoking the fast chemistry limit. 
Emphasis is placed on the ignition time and location as a function of vortex Reynolds number and initial 
temperature differences of the reacting species. In the next section we state the problem. Section 3 presents 
an analysis of the ignition regime using a combination of large activation energy (and hence Zeldovich 
number) asymptotics and numerics. Section 4 presents a simple model problem for the diffusion flame 
regime describing the viscous core region about the flame sheet, and identifies certain physical effects by 
employing asymptotic expansions and numerical solutions. The full equations are then solved numerically 
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in section 5. Comparisons between asymptotics and numerics are given through out the study. Finally, our 
conclusions are presented in section 6. 


2. EQUATION FORMULATION 


In this section the problem for the time evolution of initially unmixed species occupying adjacent 
half-planes which are then allowed to mix and react in the presence of a vortex is stated. The nondimen- 
sional equations, boundary conditions and initial conditions governing this field in cylindrical coordinates 
are 


o 

II 

V = ^-[l - exp(-r 2 /4/Sc)], 
r 

(2.1a) 


T, + jTq = V 2 T + pa 

(2.1b) 


+ = V 2 Fj - a 7 = 1.2. 

(2.1c) 


O = DaF l F 2 e~ Z£lT , 


(2. Id) 


T = F X = 1,F 2 =0 


at t = 0, r > 0, 0 < 0 < 7T, and t > 0, r 0 < 0 < n. 


(2.2a) 


T = p 7 , Fi = 0, F 2 = “7 r — <(> 1 at t = 0, r > 0, n < 0 < 2% , 

t i >ot> 


and t > 0, r % < 0 < 2 rc, 


(2.2b) 


where (U,V) are the radial and angular velocity components, respectively; V 2 is the two-dimensional Lapla- 
cian operator in cylindrical coordinates; T is the temperature; and F x and F 2 the mass fractions. The chemi- 
cal model is a one-step irreversible Arrhenius reaction, and the vortex model is an Oseen vortex with circu- 
lation T. The nondimensional parameters appearing above are the vortex Reynolds number R , defined as 
T / 2tcv where v is the kinematic viscosity; the Schmidt number Sc = vID assumed equal for both species 
and D is the species diffusion coefficient; the Zeldovich number Ze = E I (R° TJ) with E the dimensional 
activation energy and R° the universal gas constant; the Damkohler number Da defined as the ratio of the 
characteristic diffusion time scale to the characteristic chemical time scale; P the heat release per unit mass 
of F x \ the equivalence ratio <)> defined as the ratio of the mass fraction F x to the mass fraction F 2 divided 
by the ratio of their molecular weights times their stoichiometric coefficients; and finally p 7 the ratio of the 
initial temperatures. If <f> = 1, the mixture is said to be stoichiometric, if <j> > 1 it is F x rich, and if <j> < 1, it 
is F i lean. Also, if p 7 is less than one, F 2 is relatively cold compared to F u and if p 7 is greater than one it 
is relatively hot. The temperature and mass fractions were nondimensionalized by the values T and F l ooy 
respectively. The velocities were nondimensionalized by U 0 , some characteristic speed. Lengths and times 
are referred to the relevant diffusion characteristic scales l d = \ipC p U Q and l d lU 0 of the flow, 
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respeclively, where p is the density assumed constant, X is the thermal conductivity, and C p is the specific 
heat at constant pressure. Finally, the assumption of unit Lewis number is made, so that Sc = Pr, where Pr 
is the Prandtl number. 

3. IGNITION REGIME 

At time t = 0. the reaction rate is exactly zero owing to the product F ^2 = 0. For t > 0, the fuel and 
oxidizer begin to mix by diffusion, as well as by convection due to the presence of the vortex, and the reac- 
tion rate is no longer zero. For small time, it can reasonably be assumed that the effect of the reaction is 
small, and in the absence of reaction term the system (2.1) reduces to 

T, + ^-4r~ [1 - exp(-r 2 /4 1 Sc)]T 6 = V 2 J , (3.1a) 

r l 


Fj , + [1 - exp(-r 2 /4/ Sc)]F y , 0 = V 2 F ; , j = 1,2. (3.1b) 

r 1 

m 

The above system must be solved numerically subject to the conditions (2.2), and its solutions are referred I 

to as the inert solutions, denoted by f 

| 

7) = 7j(r ,9,0, Fjj s Fjj(r ,0,/), j - 1,2. (3.2) 

^ time increases, more of the combustible mixes until, at some finite time, a thermal explosion occurs | 

characterized by significant departure from the inert. I 

To analyze the ignition process, we determine the effect of the growing reaction rate by expanding I 

about the inert solution as 


T = Tj +Ze~ l T u Fj =F JtJ +Ze~ l F jtl j = 1,2, (3.3) 

and take the asymptotic limit Ze — » Substituting into (2.1), we see that the left hand sides are 0(Ze ) 
while the right hand sides are 0(Da e le T *), where T+ = max(Tj). In order to balance these terms, Da 
must therefore have the form 


Da = 


1 

pZe 


Ze iT m 

e 


Substituting into (2.1b) yields the following equation for Tj 


(3.4) 


r u + ^f £ [l-exp(-r 2 /4rS C )]r lie = V 2 7, 


7*! + Ze (T) - T. ) 

T. (7/ + Ze~ l T i) J’ 


+ [Fj,/ +Ze“ 1 F lil ] [E 2i/ +Ze 1 F 2 , 1 ]exp 


(3.5) 
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which must be solved subject to the initial and boundary conditions 

T x = 0 at t =0 y r > 0,0 <Q <2n t and t > 0, r -» «> f 0 < 0 < 2 n. (3.6) 

Once Ti is known, the mass fraction perturbations can be found directly from (2.1c). Note that the right 

hand side of (3.5) is exponentially small except where \Tj - T* \ ~ 0(Ze~ l ). There are now two cases to 

consider, depending upon the magnitude of the parameter 1 1 - p r |. In this section we only consider the case 

of nearly equal initial temperatures (P r = 1), and will postpone the. analysis for 0(1) initial temperature 
differences as a topic for future research. 

As mentioned above, only the important case of ignition for nearly equal initial temperatures of the 
fuel and oxidizer is considered. To this end, we set 


Pr»l +Ze- l pr, (3.7) 

with j3 r fixed and (7(1). With this choice, the inert temperature T, is given by 

7} = 1 + Ze~ 1 T lti (3.8) 


where T t l satisfies (3.1a) subject to the initial and boundary conditions 

(3.9a) 
(3.9b) 


T, j = 0 at / = 0, r > 0, 0 < 0 < n, and t > 0, r -» 0 < 0 < 7t, 

Tj i = B r at t = 0, r > 0, n < 0 < 2 7t, and t > 0, r -> 7t < 0 < 2 7t. 


Substituting into (3.5), together with the asymptotic choice 7. = 1, yields the following equation for the dis- 
turbance temperature 


T u + ^[1 -exp(-r 2 /4r5c)]T 1 , e = V 2 7, + F u F 2>/ e 


r, +r, 


M 


(3.10) 


subject to the conditions (3.6). This equation was solved numerically for a range of R and with 
Sc = <j> = 1. Implementation of boundary conditions are facilitated if the system is recast in Cartesian coor- 
dinates. The solution technique is a 2nd-order finite difference scheme on a nonuniform mesh. To resolve 
the structure in the core region of the field, a coordinate stretching of the form 


y P = 


C2-I 

c 2 -yc 2 


r 

ym^y c 


(3.11) 


is used, where y p e [->w, ymJ is the physical coordinate, y c e [-1,1] is the computational coordinate, 
and Ci and C 2 are adjustable constants. The smaller the quantity (C 2 - 1), the stronger the stretching. In 
the y direction, C x is either 4 or 2, and C 2 is 2. A similar stretching is used in the x direction, with C\ = 4 
and C 2 = 2. To avoid the singularity at the origin, no mesh points are placed there. The outer boundaries 
are set at 50 or 200 for and 20 or 50 for y m4X . Grid resolution studies which at least doubled the com- 
putational mesh were carried out to ensure that structures were well resolved. The resolutions required 
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ranged from a 64 2 mesh to a 256 2 mesh for large vortex Reynolds number (R = 5000). The time-stepping 
scheme is a four-stage Runge-Kutta which is formally 2nd-order but has an extended stability region making 
it accurate and robust for moderately stiff problems. All runs were performed interactively on a Cray YMP. 

As t increases, the solution for T j becomes unbounded at some finite time (t ig ) and location (x ig , y lg ). 
This characterizes the ignition regime. The special case R = 0 t in which two initially unmixed species are 
allowed to diffuse without the mixing generated by the vortex, corresponds to the results of Linan and 
Crespo (1976). To verify the solution technique described above, we give a comparison of our results to 
that of Linan and Crespo in Table I. For this case ignition takes place along a line parallel to the x-axis 
which is located at y = 0 for j$ r = 0 and resides in the hotter region for * 0. 

The effects of the vortex can be seen by examining Figures 1 and 2 and Table II. Figure 1 is a plot 
of the ignition times t lg versus the vortex Reynolds number R for three values of pr with Sc = <)> = 1. 
Table II contains selected ignition times as a function of R for the three values of of Figure 1. The 
effect of increasing from zero (i.e., p r > 1) is to enhance ignition, while decreasing from zero (i.e., 
< 1) has the opposite effect. One can see from Figure 1 and Table II that the vortex Reynolds number 
has little effect on the overall ignition time. However, ignition now occurs at a point rather than along an 
entire line. This effect is shown in Figure 2 where the ignition locations y ig versus x ig are plotted for vari- 
ous values of R and $ T = -2. Note that as R increases, the ignition location spirals clockwise towards the 
viscous core center. Ignition is seen to occur within the core for vortex Reynolds numbers/? > 70. This is 
in contrast to the case of equal initial temperatures (£ r = 0), where ignition always occurs at the origin for 
any value of R . The analogous plot for = 2 has essentially the same characteristics as Figure 2. In sum- 
mary, ignition occurs in the region of the initially hotter reactant for small R , and spirals clockwise towards 
the viscous core center as R increases. 

Figure 2 indicates that the ignition location is at its maximum from the center of the vortex at approx- 
imately R = 28. Surface plots of T } and F 1;/ shown in Figures 3a, b, respectively, display the skewed struc- 
ture of the reaction core at these conditions. The development of the hot spot for = - 2 as a function of 
increasing R is shown in the contour plots of T x and F\ j in Figures 4 and 5, respectively. The structure of 
the reaction core develops from an asymmetrically skewed center (Figures 4a-c) to an axisymmetric center 
(Figure 4d). However, the local structure of the hot spot within the reaction core is axisymmetric for any 
value of the vortex Reynolds number greater than zero. This fact is more clearly seen in Figure 6 where we 
plot T i versus y at the ignition location and time t ig corresponding to the same conditions given in Fig- 
ure 4. From this figure it is clear that once the hot spot develops within the reaction core, it is axisym- 
metric. 

For completeness, a contrast of the effects of a potential vortex distribution (V =R Sc Ir) versus the 
Oseen vortex, given in (2.1a), on the ignition locations are made. Figure 7 is a plot of y ig versus*^ for the 
potential vortex for the same conditions as in Figure 2. Note that the collapse of the spiral into the center 
occurs for R > 7, as opposed to the value R >10 for the Oseen vortex. 

The numerical solutions presented above suggests that near thermal runaway, a hot spot develops at 
(r ig ,0 t g) within the reaction core and that the structure of this hot spot is axisymmetric (see Figure 6). Thus, 
convection is not important and the local structure of the hot spot is diffusion controlled. In this case, Dold 
(1985, 1989) has shown that the proper scalings characterizing the local structure of the hot spot are given 
by 


ini hi 1 1 urn i i ii hi i j i mi iiiii i i mi mi 
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* = r ig cos 0 lft + H(t) i cos 0, y = r ig sin 0, ? + n(x) j sin 0, t = t ig - I, (3.12) 


where 


|i 2 (x) = T [ a - In (Q ig t)] = t£, (3.13) 

as t — > 0, with Cl tg - F X I F 2 t i e /1 evaluated at (r,0,f) = (r^,0 ig ,/^), and a is a constant determined by 
matching with the initial conditions and is a function of the parameters (/?,Sc,<|>). With these scalings, and 
ignoring all dependence on 5, the asymptotic form of T x in the hot spot is given by 


T } = -ln(Q* t) - ln(l + s 2 /4) - 


sjj_ 4 
1 + j 2 /4 


+ 


% 


1 +n 
2 


s 2 i 4 
1 + s 2 f 4 


In (1 + s 2 ! 4) 


+ 0 ( 



), 


(3.14) 


with n = 0 if r lg * 0 or n = 1 if = 0. Again note that the local structure of the hot spot, once it 
develops, is independent of the influence of the vortex and so is diffusion controlled. We are currently 
investigating the future time development of the hot spot, and the subsequent flame development. 


4. DIFFUSION FLAME REGIME 


Typically, a diffusion flame is characterized by a chemical reaction time that is much smaller than a 
characteristic diffusion time. Chemical reactions then occur in a narrow zone between the fuel and the oxi- 
dizer, where the concentrations of both reactants are very small. Mathematically, the assumption of very fast 
chemical reaction rates leads to the limit of infinite Damkohler number which reduces the diffusion flame to 
a flame sheet (i.e., local chemical equilibrium). This assumption significantly reduces the complexity of the 
problem since it eliminates the analysis associated with the chemical kinetics. For many flows, the assump- 
tion of local chemical equilibrium adequately predicts the location and the shape of the diffusion flame 
(Buckmaster and Ludford, 1982; Williams, 1985). Since the flame sheet model is amenable to analysis, we 
give its structure below. For finite values of the Damkohler number, equations (2.1) must be solved numeri- 
cally and this is done in the next section. 

We begin the analysis of the diffusion flame regime by defining the following conserved variables 
(see, e.g.; Williams, 1985): 


T + p/ 7 ! - p r + (1 - P/- + P)Z, 


(4.1) 


T + p/ 7 2 -P7’-i-P<j> 1 +(l — pr-|}<|> l )Z t 

_ F i - F 2 + ( jT 1 
i + <jr’ 


(4.2) 


(4.3) 
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where Z is the mixture mass fraction and satisfies the convection-diffusion equation 

Z, + ^[1 -exp(-r 2 /4rSc)]Z 6 = V 2 Z, (4.4) 

r 

silbject to the initial and boundary conditions 

Z = 1 at t = 0, r > 0, 0 < 6 < jc, and t > 0, r -» «>, 0 < 0 < Jt, (4.5a) 

Z = 0 at t ~ 0, r > 0, jr < 9 < 2n, and t > 0, r — > n < 0 < 2n. (4.5b) 

In the limit of infinite Damkohler number the flame sheet solution is given by 

Ft = 1 - (1 + <]) _1 )(1 - Z), F 2 = 0, (4.6a) 

r = Pr + p<)f 1 + (l-p 7 --p<tf l )Z, (4.6b) 

valid for Z > Z f , and 

Ft = 0. F 2 = <t> _1 - (1 + <t> _I )Z. (4.7a) 



Note that the flame location is independent of |V and p. Once Z = Z(r,0,/) is known, then the other vari- 
ables (T ,F t.F^ can be found from (4.6)-(4.7). 

As discussed in the Introduction, equation (4.4) for Z has been solved numerically by Laverdant and 
Candel (1988) and Rehm, Baum, Lozier and Aronson (1989). Our goal is not the numerical solution of 
(4.4). Rather, we propose in the following subsections a simple model problem describing the viscous core 
center in a region about the flame sheet, with the goal of identifying certain physical effects by employing 
asymptotic expansions. 


Illllili 
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4.L LOCAL ANALYSIS ABOUT FLAME SHEET 


In this section we formulate a model problem describing the effect of the vortex in a neighborhood 
about the flame sheet. Assume at time i = 0 a flame sheet exists at >y. In the absence of a vortex, the mix- 
ture mass fraction evolves according to (4.4) with R = 0, and the analytical solution can be written in terms 
of an error function. At any instantaneous time t Q , the solution across the flame sheet is easily seen to be a 
linear shear profile (Figure 8), given by 

Z(.r,Q,t 0 ) = + Z/(0(y - y f (t 0 )l (4.9) 

which is only the first-order Taylor series expansion of the mixture mass fraction profile Z about the flame 
sheet location found implicitly by (4.8a). Here, the prime denotes differentiation with respect to y. The 
model problem assumes that the time scale associated with the vortex is faster than that of diffusion, so that 
on this time scale the linear profile (4.9) can be considered steady. In addition, we assume that the flame 
sheet resides close to the origin. Thus the following scales can be introduced 

t = t 0 + £ f , y=Vey, r=Ver, <|>=1 + V£$, Z = y + Viz, (4.10) 

with £ << 1, and the vortex is turned on at time / = 0. With these scalings the linear shear profile, in terms 
of Z , becomes 


Z{r,Q,0) = -^+Z;« o )[? ^k l+ k 2 y. (4.11) 

Note that k l = 

Z f = 1 / 2 at t 
the flame sheet 
Z becomes 


-<j)/4 - Zf(t 0 )y f {t 0 ) measures the instantaneous deviation of Z from its stoichiometric value 
= 0, and k 2 = Zf(t 0 ) measures the instantaneous gradient of the mixture mass fraction across 
. The convection-diffusion equation (4.4) in terms of the perturbation mixture mass fraction 


Z ; + [1 - exp(-r 2 /4f Sc)]Z e = % 2 Z, (4.12) 

where V 2 is the Laplacian operator in the transformed plane. The perturbation equation (4.12) is to be 
solved subject to the initial and boundary condition given by the linear shear profile (4.11). 

To examine the time evolution of the flame sheet about the origin when the vortex is present, we 

define 


Z(r,0,O = ifc 1 + k 2 [f sin© + Z(r, 0 , 0 ], (4.13) 

where Z is the disturbed mass fraction due to the presence of the vortex. Upon substituting (4.13) into 
(4.12), we find that Z satisfies the convection-diffusion equation 
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Z- + [1 - exp(-r 2 /4/ Sc)] (f cosO + Z 9 ) = ^ 2 Z, (4.14) 

f 


subject to homogeneous initial and boundary conditions. Note that (4.14) is non-homogeneous, and hence 
does not have a similarity solution. Here, we introduce the variable 


f 



(4.15) 


which is the similarity variable in the absence of the vortex, so that Z ■ = Z(s ,0,r). In the transformed plane 
Z satisfies 

t Z- - jZ s + - y- [1 - exp(-s 2 /45c)](s VT cos9 + Z 0 ) = P 2 Z. (4.16) 


This equation is separable in space and time and hence has the simple analytical solution 

Z = V7 [Z e (s)cos0 + Z o (s)sin0], (4.17) 


where 

-s±')+ ^| £ [l-exp(-s 2 /45c)](Z 0 + s) = z"+-Z € - , (4.18) 

2 s S 

\(Z 0 -sZ 0 )- ^[1 -exp(-i 2 /4Sc)]Z e =2"+-Z^- \z o , (4.19) 

2 s l s s 

subject to the boundary conditions 

z e (0) =.Z 0 (0) = Z e (oo) = X (oo) = 0. (4.20) 

Here, primes denote ordinary differentiation with respect to the similarity variable s , Z e is the amplitude of 
the even mode, and Z a is the amplitude of the odd mode. Note that the above solution is valid for all times 
and vortex Reynolds numbers. 

The system (4.18-4.19) subject to (4.20) was solved numerically using a 2nd-order finite difference 
scheme with an appropriate stretching in s. The solution to this system for a range of vortex Reynolds 
numbers is given in Figure 9. As the vortex Reynolds number increases, the number of oscillations in each 
of the components Z { and Z 0 also increases while moving away from the origin, thus establishing a core 
region. The composite solution (4.13) is plotted in Figure 10 for t = 1, 0 = ±n/ 2, kj - 1, =0, and 

Sc = 1. From this figure the extent of the core region is seen to grow as R increases, and this core, once 
established, grows like V7 (see (4.17)) as noted previously by Maible (1985). 
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4.2. LARGE TIME BEHAVIOUR 


To investigate the local structure of the model problem described in the previous section for large 
times, let t —> <», so that s — > 0 for fixed r. Hence from (4.17) we have 

Z = lim VT [Z,(r)cos0 + Z 0 (s)sinG) 


= lim — [Z,(j)cos0 + Z„(s)sin0] 

s — *0 S 


lim — [Z e (s)x + Z 0 (s)y ] 

s — >0 S 


= a x x + b x y , 


(4.21) 


where from (4.18) and (4.19) it is easy to show that the asymptotic behaviour as s — > 0 is 

Z, = a x s + 0(s 3 ), Z 0 = b x s + 0(s 3 ). (4.22) 

Here, a x and b x are to be found numerically. Substituting (4.22) into (4.13) yields 

Z=k x +k 2 [a x x + (l+b x )y]. (4.23) 


Comparing (4.23) to (4.11), the effect of the vortex is seen to establish a new equilibrium mixture mass 
fraction profile within the core region. 

Another quantity of interest is the magnitude of the gradient of the mixture mass fraction, defined as 

) 1/2 


G = 


f az 

2 

az ' 

2 ' 

dx 

v J 

+ 

l*J 

J 


(4.24) 


Note that G is a function of R , § and Sc . Substituting (4.23) into (4.24), we see that the magnitude of the 
gradient in the center of the core is given by 

G =1*2 1 [a? + (l + *i) 2 ] t2 . (4.25) 

Figure 11 is a plot of G versus R for Sc = <\> = k 2 = 1. Note that as /?-><», G 0 which implies that 
there are no gradients within the center of the core, and hence no diffusion across it. Thus, the core is com- 
pletely established. From this figure it is clear that the core is 95% established at R ~ 50. This is in agree- 
ment with Figure 10 where a well defined core is clearly visible at R = 100. 
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4.3. LARGE VORTEX REYNOLDS NUMBER ASYMPTOTICS 


To investigate the local structure of the model problem described in Section 4.1 at large vortex Rey- 
nolds numbers with R » 1, we begin with the equation 

Z i + ±l\- exp (-R r i /4t)]Z e = $ 2 Z, (4.26) 

which is (4.12) with r scaled as r = 7 SSc R . The above equation is solved subject to the initial and boun- 
dary conditions 

Z =k\+ k 2 y at t > 0, F -> °o, 0 < 0 < 2 7t, and t = 0, r > 0. 0 < 6 < 2n. (4.27) 


OUTER INVISCID SOLUTION - Z° 


In the limit R -» «>, an outer inviscid region exists and is governed by the following inviscid version of 
(4.26). 


Z?+±Z$= 0, (4.28) 

r 

where Z° denotes the mixture mass fraction in the outer inviscid region. This equation is valid for 
r = 0(1) and t <.0{R), so that the effect of a potential line vortex or an Oseen vortex is the same. The 
solution is easily found by the method of characteristics and is given by 


Z° - k } + k 2 r sin0 o , 


(4.29) 


where 0 O is the characteristic 


0 O = 0- 


t 



(4.30) 


INNER SOLUTION - Z l 


As 7 — > 0, the solution (4.29) becomes singular so the viscous terms in (4.26) must be retained. We intro- 
duce the similarity type variable 
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(4.31) 

where 5 << 1 and will be chosen in the course of the analysis. Defining Z 7 to be the inner variable, and 
expanding 



Z ] = k x +k 2 'TsT Z(s,Q), 


(4.32) 


we see that Z satisfies the convection-diffusion equation 

l A(2 - s2 s )+ -ij [1 - exp (-/? S* 2 /4)]Z 9 = — 1 ~— V 2 Z. (4.33) 

8 s l o Sc R 


To facilitate matching with the outer inviscid solution (4.29), we change to a new (non-orthogonal) coordi- 
nate system in which the characteristic 0 O is one of the coordinates. Thus, the above equation in terms of s 
and 0 O becomes 


l A(Z - s2 s ) - — ^-rexp(-/? 8 j 2 /4)Z 0 


Sr 2 
1 


85c R 


A 2 3 

3r 8 s 3 99o 


A 2 3 

ds + 8s 3 30o 


1 

+ s 2 30o I (Z ' 


(4.34) 


To balance both sides we must choose 


1 J_ 

8 RSc S 2 ’ 


(4.35) 


or, in terms of 8, 


8 = (RScT 113 . (4.36) 

For this scaling, exp( - R Si 2 / 4) is exponentially small in the limit /?—>», and so again the effect of a 
potential line vortex or an Oseen vortex is the same. The leading order equation in a 8-expansion is given 
by 


'A(2-sZ s ) = ±-2 %% . (4.37) 

This equation can be solved by separation of variables and, after matching with the outer solution (4.29), the 
solution is found to be 


-4 

3s 6 


2 = s exp 


sin0 o . 


(4.38) 
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The inner viscous solution is now given by 


Z = k\ + kz s exp 



sin 9 0 . 


(4,39) 


COMPOSITE SOLUTION 


The composite solution 


Z = ki + kz r exp 


-4(5Q 3 

37 * 


sin 




(4.40) 


is valid for times t << 5“ l . Figure 12 displays the composite solution Z versus j for 8 = ±n/2, with 
ki = 0.5, k 2 = 1, 5 = 0.05, and t = 1. Note the agreement with the numerical solution given in Figure 10. 

The magnitude of the gradient of the mixture mass fraction can be obtained by substituting (4.40) into 
(4.24). Then as 5 -» 0, the leading order result yields 


G = 





(4.41) 


Figure 13 is a plot of G versus j for 8 = 7t/2, k 2 = 1, and 8 = 1/ 10. Note that the maximum amplification 
is about 5~'. The limit R -» 0 corresponds to G -» 0, and so the large spikes have been created by the vor- 
tex. . 

The core radius growth can be determined from (4.31) and (4.36), and is found to be 


r 

V(/?e Sc)- m t 


= constant t 


(4.42) 


or, in terms of dimensional quantities (denoted by stars), 


Vr^'V'Y 


= constant. 


(4.43) 


This is the main result of Marble (1985), who showed constant = 0.5092 + 0(VD* / T* ). Finally, we note 
from (4.43) that the spreading of the viscous core due to the vortex field is (T ID )^ 3 times as large as 
that obtained from diffusion alone in the absence of the vortex, as previously pointed out by Marble (1985) 
and Karagozian and Marble (1986). 
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5. NUMERICAL SOLUTION 


In this section we present selected numerical results to the full system (2.1) subject to (2.2). This sys- 
tem was solved by a standard finite difference scheme, as previously described in Section 3. To illustrate the 
numerical solution of (2. 1-2.2) of the continuous evolution from nearly frozen flow to near equilibrium flow, 
we produce here one result corresponding to the case Sc = <j> = 1, R = 28, Ze - 30, Da = e 30 / 30, (3=1, 
and p r = 1 - 2/30. This case corresponds to the asymptotic results of Section 3 as displayed in Figures 3- 
6. Figure 14a is a plot of the time slices of the temperature profile 7\ while Figure 14b is a plot of the time 
slices of the reaction rate term ft. Time increases from the bottom left comer to the top right comer. In 
Figure 14a, a hot spot develops initially within the viscous core (frame 1) and rapidly develops into an iso- 
lated, almost circular flame which grows as time increases (frames 2-8). At a later time (frame 7) two 
diffusion flames are clearly visible at the edges of the plot, and as time increases further (frame 8) they 
move towards the viscous core. Eventually the two diffusion flames will merge with the expanding, almost 
circular flame located in the viscous core region. Figure 14b shows corresponding results for the reaction 
rate term. Initially (frame 1) a single isolated point in the rate term is visible. As time increases (frames 2- 
8) this single point develops into an almost circular ring and spreads in time, leaving within it the burnt 
core. The reaction rates of the two diffusion flames are given in frames (7-8). 


6. CONCLUSIONS 


The distortion of flames in flows with vortical motion has been examined by means of asymptotic 
analysis and numerical simulation. The model consisted of a constant-density, one-step, irreversible 
Arrhenius reaction between initially unmixed species occupying adjacent half-planes which were then 
allowed to mix and react in the presence of a vortex. The continuous evolution of the temperature and mass 
fraction fields from initially unmixed to near equilibrium flow was followed. Emphasis was placed on the 
ignition time and location as a function of vortex Reynolds number and initial temperature differences of 
the reacting species. In the ignition regime, the case of near equal initial temperatures was considered. The 
effect of increasing the initial temperature ratio p r slightly from unity was to enhance ignition for any fixed 
vortex Reynolds number R , while decreasing (3 r slightly from unity had the opposite effect. However, for 
fixed Pr, increasing the vortex Reynolds number from zero had little effect on the overall ignition time. The 
ignition location occurred in the region of the initially hotter reactant for small R , and was seen to spiral 
clockwise towards the viscous core center as R increases. Finally, numerical solutions of the ignition equa- 
tions indicated that the hot spot was both axisymmetric and diffusion controlled for any non-zero vortex 
Reynolds number. Thus, the asymptotic analysis of Dold (1985, 1989) could be used in describing this local 
structure. The case of different initial temperatures, and the future time development of the hot spot, will be 
considered in a future manuscript. 

In the diffusion flame regime, a simple model problem was proposed to investigate the effect of the 
vortex on the flame sheet. The model consisted of assuming that the time scale of the vortex was faster than 
that of diffusion alone, so that on this time scale appropriate asymptotic expansions were used in describing 
the interaction about the flame sheet. This model problem allowed for the full vortex flow field without 
making the strain-shear approximation of Marble (1985). Many of the conclusions of Marble (1985) were 
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verified, including the similarity rule for the core radius growth in the limit of large vortex Reynolds 
number. 
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Figure I. Plot of lire ignition time 1^ versus vortex Reynolds number R for Sc = $ = 1 and (a) j^- = -2 
(b) = 0, and (c) ^ = 2. 
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Figure 2. Plot of die ignition locations versus Xj, for various values of the vortex Reynolds number 
with Sc - - 1 and - - 2. 
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Figure 3. Surface plot of (a) Tj and (b) F l ( for R = 28.0, (V[- = -2, and Sc 
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Figure 3. Surface plot of (a) Tj and (b) F\j for R! = 28.0, j^p = -2, and Sc 
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(a) 



Figure 5. Conlour ploi of Fi 4 for j^- - - 2, Sc - 4* - 1 and (a) R = 5, (b) R - 20, (c) R = 28, and (d) 
R - 100, 
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Figure 7. Plot of the ignition locations y^ versus x^ for various values of tire vortex Reynolds number 
widi Sc = <|> = 1 . jij- = - 2, and using die potential vortex defined by V = RSc/r. 



Figure 8. Schematic showing the assumed instantaneous flow field. 





Figure 10. Plot of Z versus s for t - 1, 0 - ± rt/2, k, = 0, k 2 = I, Sc = 1 and various values of the vorlex 
Reynolds number. 



R 


Figure 11. Plot of G 


versus vorlex Reynolds number for Sc = $ = k 2 = 1. 
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Figure 14. Plot of the time slices of (a) the temperature profile T and (b) the reaction rate term Cl. Time 
increases from the bottom left corner to top right comer. 
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